Nucleon-Nucleon Phase Shifts and Pairing in Neutron Matter and Nuclear Matter 
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We consider ^So pairing in infinite neutron matter and nuclear matter and show that in the 
lowest order approximation, where the pairing interaction is taken to be the bare nucleon-nucleon 
(NN) interaction in the ^So channel, the pairing interaction and the energy gap can be determined 
directly from the ^Sq phase shifts. This is due to the almost separable character of the nucleon- 
nucleon interaction in this partial wave. Since the most recent NN interactions are charge-dependent, 
we have to solve coupled gap equations for proton-proton, neutron- neutron, and neutron-proton 
pairing in nuclear matter. The results are, however, found to be close to those obtained with 
charge-independent potentials. 
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Recently, there has been renewed interest in the pair- 
ing problem in neutron matter and neutron-rich nuclei. 
The superfluid properties of neutron matter is of impor- 
tance in the study of neutron stars [0, while pairing in 
neutron-rich systems is of relevance for the study of heavy 
nuclei close to the drip line 1^ and the light halo nuclei 
Much effort has gone into calculating the superfluid 
energy gap in dilute neutron matter . Most of these 
studies, e.g., those of Refs. have been carried out 

using pairing matrix elements given by the bare nucleon- 
nucleon (NN) interaction. Many of the same authors 
have calculated the ^5*0 gap in nuclear matter, which has 
also been the subject of recent relativistic formulations 
of the pairing problem |p|-pl| . 

Even though it is a long time since Clark et al. fl^ j 
showed that the effects of density and spin-density fluc- 
tuations must be included in the pairing interaction, and 
there has been much progress in that direction recently 
[^,0 , we will here focus on the situation at the level of 
the bare interaction. In this lowest-order approximation 
to the problem it has been found that results for the ^S'o 
energy gap in neutron matter and in nuclear matter are 
almost independent of the choice of NN interaction. We 
aim at explaining how this can be understood directly 
from the measured properties of the free NN interaction. 
Moreover, although a relation between the pairing gap 
and NN phase shifts was obtained almost forty year s ag o 
by Emery and Sessler jl5| (see also Hoffberg et al. [[l6|), 
in this work we wish to focus on the near interaction in- 
dependence of the results for the energy gap at the Fermi 
level, and try to explain this from the NN scattering data 
directly. 

The energy gap in infinite matter is obtained by solving 
the BCS equation for the gap function A(fc). 
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where V{k, k') is the bare momentum-space NN inter- 
action in the ^S'o channel, and E{k) is the quasiparti- 
cle energy given by E{k) = \J{eik) — e{kF)Y + A(fc)2, 



where e(fc) is the single-particle energy of a neutron 
with momentum fc, and kp is the Fermi momentum. 
Medium effects should be included in e(fc), but we will 
use free single-particle energies e(fc) = fc^/2m, where m 
is the neutron rest mass, to avoid unnecessary compli- 
cations. And in neutron matter, at least at the densi- 
ties considered here, Brueckner-type calculations [Q in- 
dicate that in-medium single-particle energies do not dif- 
fer much from the free ones. The energy gap is defined 
as Ai? = A(fci?). Eq. (|l|) can be solved by various tech- 
niques, some of which are described in Refs. In Fig. 
|l| we show the results for A p obtained with the CD-Bonn 
potential (full line) ||l^, the Nijmegen I and Nijmegcn II 
potentials (long-dashed line and short-dashed line, re- 
spectively) [HI . The results are virtually identical, with 
the maximum value of the gap varying from 2.98 MeV for 
the Nijmegen I potential to 3.05 MeV for the Nijmegen 
II potential. The same insensitivity of the energy gap 
with respect to the choice of NN interaction was found 
in Refs. We will now discuss how these results can 
be understood from the properties of the NN interaction 
in the ^5*0 channel. 

A characteristic feature of ^5*0 NN scattering is the 
large, negative scattering length, indicating the presence 
of a nearly bound state at zero scattering energy. Near 
a bound state, where the NN T-matrix has a pole, it can 
be written in separable form, and this implies that the 
NN interaction itself to a good approximation is rank- 
one separable near this pole |l^. Thus at low energies 
we can write 



V{k,k') = \v{k)v{k'), 



(2) 



where A is a constant. Then it is easily seen from Eq. (|l|) 
that the gap function can be written as Apv{k), where 
Ai? is the energy gap. Inserting this form of A(A;) into 
Eq. (|l|) one obtains 
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which shows that the energy gap Ap is determined by 
the diagonal elements Xv^{k) of the NN interaction. The 
crucial point is that in scattering theory it can be shown 
that the inverse scattering problem, that is, the determi- 
nation of a two-particle potential from the knowledge of 
the phase shifts at all energies, is exactly, and uniquely, 
solvable for rank-one separable potentials ||l^,|2^. Fol- 
lowing the notation of Ref. Il^ we have 



+ k\ sin(5(fc) 



-a(k) 
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for an attractive potential with a bound state at energy 
E = —k\. In our case kb = 0. Here 5{k) is the ^Sq 
phase shift as a function of momentum k, while a{k) is 
given by a principle value integral: 



a{k) = -P / dk 

I" J-oo 



-fc' 



(5) 



where the phase shifts are extended to negative momenta 
through S{—k) = —S{k). Eqs. and (0) can also be 
rewritten in terms of the Jost function [ pO| as done in 
Ref. 111. 

From this discussion we see that Aw^(fc), and there- 
fore also the energy gap Af, is completely determined 
by the ^Sq phase shifts. However, there are two obvious 
limitations on the practical validity of this statement. 
First of all, the separable approximation can only be ex- 
pected to be good at low energies, near the pole in the 
T-matrix. Secondly, we see from Eq. that knowledge 
of the phase shifts S{k) at all energies is required. This 
is, of course, impossible, and most phase shift analyses 
stop at a laboratory energy Ei^h — 350 MeV. Strictly 
speaking, the rank-one separable approximation to the 
^5*0 interaction breaks down already where the phase 
shift changes sign from positive to negative at i?iab ~ 248 
MeV, corresponding to a single-particle momentum of 
k « 1.73 fm~^. However, at low values of kp, knowledge 
of v(k) up to this value of k may actually be enough to 
determine the value of Ap, as the integrand in Eq. (||) is 
strongly peaked around kp. We therefore found it worth- 
while to try to calculate the energv gap directly from the 
^5*0 phase shifts using Eqs. (||)-(||). A possible improve- 
ment to the rank-one separable approach for potentials 
which change sign is discussed by Kwong and Kohler |l] . 

The input in our calculation is the ^5*0 phase shifts 
taken from the recent Nijmegen phase shift analysis [ p2| . 
We then evaluated Xv'^(k) from Eqs. (||) and (||), using 
methods described in Ref. |2^ to evaluate the principle 
value integral in Eq. (H). Finally, we evaluated the energy 
gap Ap for various values of kp hy solving Eq. (||). Nu- 
merically the integral on the right-hand side of this equa- 
tion depended very weakly on the momentum structure 
of A(fc), so in our calculations we could take A(fc) ^ Ap 
in Eq. (^, and thus it became an algebraic equation for 
the energy gap A p . The resulting energy gap is plotted in 
Fig. 1^ (dashed hue) together with the gap obtained with 
the CD-Bonn potential (full line). As the reader can see, 



the agreement between the direct calculation from the 
phase shifts and the CD-Bonn calculation oi Ap is very 
good, even at densities as high as kp = 1.4 fm^^. The 
energy gap is to a great extent determined by the avail- 
able ^Sq phase shifts. In the same figure we also report 
the results (dot-dashed line) obtained using the effective 
range approximation to the phase shifts: 



k cot S{k) = -— + -^rofc^. 



(6) 



where oq = —18.8 ± 0.3 fm and tq = 2.75 ± 0.11 fm are 
the singlet neutron-neutron scattering length and effec- 
tive range, respectively. In this case an analytic expres- 
sion can be obtained for Aw^(fc), as shown in Ref. pOj]: 
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and 



-2/aoro, and (3i : 

are the two roots of the quadratic equa- 



(8) 



The phase shifts using this approximation are positive at 
all energies, and this is reflected in Eq. (^ where Xv^ {k) 
is attractive for all k. From Fig. ^ we see that below 
kp = 0.5 fm"^ the energy gap can with reasonable accu- 
racy be calculated with the interaction obtained directly 
from the effective range approximation. One can there- 
fore say that at densities below kp = 0.5 fm~^, and at the 
crudest level of sophistication in many-body theory, the 
superfluid properties of neutron matter are determined 
by just two parameters, namely the free-space scattering 
length and effective range. At such densities, more com- 
plicated many-body terms are also less important. Also 
interesting is the fact that the phase shifts predict the 
position of the first zero of A(k) in momentum space, 
since we see from Eq. (^ that A{k) = Apv{k) = first 
for S{k) = 0, which occurs at i?iab ~ 248 MeV (pp scat- 
tering) corresponding to fc « 1.73 fm~^. This is in good 
agreement with the results of Khodel et al. Q . In Ref. 
Q it is also shown that this first zero of the gap function 
determines the Fermi momentum at which Aj? = 0. Our 
results therefore indicate that this Fermi momentum is 
in fact given by the energy at which the ^Sq phase shifts 
become negative. 

The calculation of the ^Sq gap in symmetric nuclear 
matter is closely related to the one for neutron matter. In 
fact, with charge-independent forces, like the older Bonn 
potentials, and free single-particle energies one would, 
of course, obtain exactly the same results. However, the 
new potentials on the market are charge-dependent, in 
order to achieve high quality fits to both np and pp scat- 
tering data, and therefore we must solve three coupled 
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gap equations for neutron-neutron (nn) , proton-proton 
(pp), and neutron-proton (np) pairing |2J]: 



(9) 



where z=nn, pp and np, and the quasiparticle energy is 
stih given by E{k) = ^(e(fc) - e(fci.))2 + A(fc)2, but the 
energy gap is now given by 

A(fc)2 = A„„(fc)2 + App(fc)2 + A„p(fc)2. (10) 

Thus the equations are coupled through their common 
energy denominator. The ^Sq pp and nn interactions are 
very nearly identical, so the set of equations above can be 
reduced to two: one for the nn (or pp) gap and one for the 
np gap. Solving these equations, both with the CD-Bonn 
potential and with the phase shift approximations we get 
the results shown in Fig. ^. For comparison we have 
in the same figure plotted the results for pure neutron 
matter with the CD-Bonn potential (dashed line). From 
the figure it is clear that the phase shift approximation 
works well also in this case, and that the gap in symmet- 
ric matter is not very different from the gap in neutron 
matter. As could be expected, the results are very close 
to those obtained earlier with charge-independent inter- 
actions 

In summary, we have shown that in infinite neutron 
and nuclear matter, owing to the near rank-one separa- 
bility of the NN interaction in the ^Sq partial wave, we 
are able to compute the ^5*0 pairing gap directly from the 
NN phase shifts. This explains also why all NN poten- 
tials which fit the scattering data result in almost identi- 
cal ^5*0 pairing gaps. This is the main result of this work. 
However, it should be mentioned that this result is not 
likely to survive in a more refined calculation, for instance 
if one includes density and spin-density fluctuations in 
the effective pairing interaction like in e.g., Refs. p3| , p^ . 
Other partial waves will then be involved, and the sim- 
ple arguments employed here will no longer apply. Our 
reasoning here applies also only to a partial wave where 
the T-matrix (almost) has a pole, and we have neglected 
the fact that the phase shifts become negative at higher 
energies. As a curiosity, we have found that at Fermi mo- 
menta below 0.5 fm~^ the pairing gap is even determined 
by two parameters only, the effective range and the scat- 
tering length. Also, we have pointed out that since the 
new NN interactions are charge dependent, one has to 
consider three coupled gap equations for ^5*0 pairing in 
nuclear matter. The final result though is very nearly 
the same has what one obtains with charge-independent 
interactions. 

We are much indebted to J. W. Clark, and E. Osnes 
for many valuable comments and discussions. 
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FIG. 3. ^So energy gap in nuclear matter calculated with the CD-Bonn potential compared with the direct calculation from 
the ^So np and pp phase shifts. Also shown are the results for neutron matter with the CD-Bonn potential. 
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